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Рассматривается задача идентификации дебита скважин нефтяного месторождения с учетом дополнительной информации о 
дебитах (приемистости) соседних скважин окружения, дополнительной априорной информации и экспертных оценок параме- 
тров модели притока жидкости, представленных непараметрическими моделями. Приводятся результаты статистического моде- 
лирования по определению точности предложенных оценок продуктивности скважин и пластового давления. 

Актуальной проблемой мониторинга и опера- ний в связи с неполнотой, неоднородностью про- 
тивного управления работой скважин месторожде- мысловых данных, результатов исследований сква- 
ний углеводородов является задача идентифика- жин, наличия различного рода ошибок, отсутстви- 
ции гидродинамических параметров (дебита неф- ем достоверной информации о моделях взаимодей- 
ти, жидкости, воды, пластового давления, фильт- ствия скважин, моделей дополнительных априор- 
рационных параметров пласта и т. д.), которая за- ных сведений и экспертных оценок, 
ключается в построении оптимальных, в смысле в этой связи актуальным является идентифика- 
заданных критериев качества, математических мо- ция гидродинамических параметров пластов с не- 
делей на основе промысловых данных, результатов пользованием метода интегрированных моделей 
исследований скважин с использованием дополни- [3—5], который позволяет учитывать различную 
тельной априорной информации, экспертных оце- неоднородную, дополнительную априорную ин- 
нок, накопленного опыта и знаний [1-3]. формацию, дает оптимальное решение задачи 

Использование классических методов иденти- идентификации с обеспечением устойчивости ре- 

фикации гидродинамических параметров [1,2] ча- шения, согласованности исходных данных и до- 
сто связано с проблемами низкой точности реше- полнительных априорных сведений. 
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В данной работе рассматривается интегриро- 
ванная система моделей дебита скважины вида: 

= = /А МО, %'(*,)) +Щ), 

■ д 3 (0 = / ѵ (<і(0) +г іА0,4 = 1 , т , (1) 

а * (1 , ) = / Зк («а- (4 )) + ѵ а (4 ) Д = 1, I, г = 1, п, 

где 4(1), Р 3 (1 ,) - измеренные значения дебита и за- 
бойного давления исследуемой скважины в момент 
времени 1{, 4,(1), ./=1 ,т - значения дебитов либо 
приемистости (для нагнетательных) скважин окру- 
жения, взаимодействующих с исследуемой сква- 
жиной; а к (1), к=\,1 - дополнительные априорные 
сведения о неизвестных значениях параметров 
а(1)=(а х (1),Ог(1),...,а к (1)) дебита исследуемой сква- 
жины к моменту времени ^ у- 1,т,/ зь /с=1,/ — 
модели дебита исследуемой скважины, дебитов 
(приемистости) скважин окружения и дополни- 
тельных априорных сведений; %(1),ц(1),ѵ(1) - слу- 
чайные величины, представляющие погрешности 
измерений притока жидкости, забойных давлений, 
ошибки, связанные с выбором модели дебита сква- 
жины и дополнительных априорных сведений. 

Следует отметить, что модель дебита исследуе- 
мой скважины / выбирается из физических зако- 
номерностей процесса нефтегазодобычи [1, 2]. 
Это, как правило, известные функции, заданные с 
точностью до параметров а(1), которые предста- 
вляют фильтрационные свойства и энергетическое 
состояние нефтяных пластов, призабойной зоны 
скважины и т. п. 

Зависимости / у , і=\,т и / ь к=\ ,1 - это неиз- 
вестные однозначные функции, априорная инфор- 
мация о которых отсутствует. Можно предполагать 
об их ограниченности, непрерывности, существо- 
вании производных и т. д. Возможно наиболее про- 
стое их представление в виде 

Л №(!>)) ~ ч(1,\ (2) 

/за(«а(4))-«а (О (3) 

в окрестности некоторых заданных значений 4(1) и 
а!(1) [5]. 

Процедура идентификации интегрированной 
системы моделей (1) с учетом приближений (2, 3) за- 
ключается в решении двух оптимизационных задач: 

а*( 1 ,И) = аг§тіп(Ф(а(0) = 

о(0 

= ^ 0 (а(^))+^ ^ (а(^))+^ 2 (а(^)), (4) 


ли дебита скважины в момент времени и упра- 
вляющих параметров 


А = (А„А 2 ,А 3 ); Л(«(*))= <7*-/і (а(і),Р 3 ) 


•/.(«(О) = Ё 


7=1 


Я,- А (а(1),Р 3 ) 


Л(«(0) = Х а ( 1 ,)~а( 1 ) 


- частные квадратичные функционалы качества 
моделей дебита скважины, дебитов (приемистости) 
скважин окружения, дополнительных априорных 


сведений 


ссі,і = 1 ,п, К , = сНа§(К 



диагональная матрица весовых функций, определяю- 
щая значимость значений дебитов скважин в различ- 
ные моменты времени истории разработки относи- 
тельно значений дебитов в текущий момент времени 1, 


ттх л- -4 ч ,/(«°(0,^(0)-?Д-), . , , 

Щ, = с 1 ш§(К( — - — )ѵѵ,( ), і = 1 ,п). 


^2 ,=(Иа8(К,(!-—^-)щ(— 

К 


(О -«;(?,.) 

К 


), І = 1 , т) 


- матрицы, в которых весовые функции щ, щ 
определяют значимость дебита скважин окруже- 
ния и дополнительных априорных сведений; а,(1), 
І=\,т - начальные приближения параметров моде- 
ли дебита скважины в момент времени I; 
И=(Н { ,к ъ Н) - управляющие параметры; 

Я=(я(іі), *=!,«)> 

/і (а (О, Р 3 ) = (/, (а (1), Р (1, ), 1 = 1, и). 


Я і = (/,(4)> і = \,п), ) = \,т 

- векторы значений дебитов жидкости исследуе- 
мой скважины (фактические и вычисленные на ос- 
нове модели/) и дебитов (приемистости) скважин 
окружения; [|А[| - норма вектора Х\ 

Весовые функции (ядра) К,щ,щ обладают сле- 
дующими свойствами [6]: 

\ѵ(х-и)/к) — > 1, к -> со; хѵ(х-и)І И) — > О, А—» 0. (6) 

Для получения оценок управляющих параме- 
тров А* (5) используется контрольный участок исто- 
рии разработки длительностью 1„_„1 п _^ ь ...,1 п с объе- 
мом выборки г. 


/?* = аг§тіп/ 0 (А) = 

к 


Анализ точности оценок продуктивности скважины 
и пластового давления 


= ^(ч(о-№У(1,мр;т 2 , (5) 

і-п-г 

где аг§тіпДх) обозначает точку минимума х функ- 
ции/; а(1,к). А* - оценки вектора параметров моде- 


Анализ точности оценок продуктивности сква- 
жины и пластового давления проводился методом 
статистического моделирования с использованием 
линейной стохастической интегрированной систе- 
мы моделей 
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Технические науки 


д\1,) = чУУУ + с Ш)) = 

= (а, У ) + а 2 • Р 3 (1. )) • (1 + (( )), 

■ д(і І ) = к-д(і і )(\ + с 2 Щ)), 

аі(г,) = а,У)(1 + с 3 ѵ 1 (?,)Х 

а 2 (1, ) = а 2 (1 + с 4 ѵ 2 (г, )), і = 1, л, 

дебита центрального ряда трехрядной схемы распо- 
ложения добывающих скважин, где коэффициент 
взаимодействия к нагнетательной скважины с при- 
емистостью д(і) и добывающей с дебитом д(і) рассчи- 
тывался дж однородного пласта при поршневом вы- 
теснении из него нефти водой по формуле [7. С. 172] 
к = (2С/В + 1)/п и «1,3 

при С/В= 1 ,42 и числе нагнетательных скважин л„=3. 

В (7) дебит скважины д(1) представлен линей- 
ной зависимостью (модель линейной фильтрации 
жидкости), в которой параметр а { (1)=Р ш (1)- а ъ где 
Р„ л (1) — пластовое давление, в течении 10 мес. раз- 
работки, изменялось по линейному закону 
(Р„ л (І)=а+Ь-і, /=1,л, а= 20, Ь= 1, л=10), а,=4 - коэф- 
фициент продуктивности скважины постоянен, 
1=і — номер месяца разработки; 4(1), і=1,п, 4(1), 
і=\,п - имитируемые значения дебитов исследуе- 
мой и приемистости нагнетательной скважины; 
а х (1), а 2 (1), /= І ,п - имитируемые значения допол- 
нительных априорных сведений о параметре а х (і) и 
коэффициенте продуктивности скважины а{, 
4,(1 •), г)(1, ),Ѵ](і), 7=1,2 - случайные величины, полу- 
ченные с использованием датчика псевдослучай- 
ных чисел, распределенных по нормальному зако- 
ну с нулевыми средними значениями и единичной 
дисперсией; с р у- 1,4- некоторые константы, пред- 
ставляющие относительный уровень ошибок. За- 
бойное давление скважины задавалось постоян- 
ным: Р 3 (і)=Р 3 = 13 МПа. 

Для модели (7) процедура определения оценок 
параметров а(1,к)=(а\(1,к),а 2 (1,к)) (4) сводится к 
решению системы линейных уравнений: 

(р' К,р + к ■ р'іѴ^р + ]Гж 2і )-а(і,к) = 

і = 1 

= (Р 1 т К,д , + Р^1Ѵ 1 д + ^1Ѵ 1І а І ), ( 8 ) 

1=1 

где Р]=(/ р /=І ,п, 7=1 ,2) - матрица размерности 
(их 2) с элементами /ц=1,/ а = Р 3 (1), /'= 1,«; К„Ж х ,\Ѵ ѵ — 
диагональные матрицы (4), (5) с весовыми функ- 
циями вида (6): К(и)=ехр(-и 2 /2); д=(д(1), і=\,гі), 
Ч=(<І(і)> *=1,й), а=(а 3 (1),]=\,2), і=\,п - векторы 
имитируемых дополнительных априорных данных 
и экспертных оценок параметров дебита скважи- 
ны; Р 3 =(1,Р 3 У - вектор-столбец. 

Следует отметить, что при значениях управляю- 
щих параметров к 2 ,к 2 =0 дополнительная априорная 
информация не учитывается, и процедура оцени- 
вания параметров (8) совпадает с взвешенным ме- 
тодом наименьших квадратов 


(Р 1 т К І Р 1 )-а(і,к 1 ) = Р^К І д\ ( 9 ) 

Для определения оптимальных значений упра- 
вляющих параметров Н (5) использовался метод 
деформированного многогранника [8]. Начальные 
приближения параметров а\і) в (8) рассчитыва- 
лись по регуляризированному методу наименьших 
квадратов [9] 

(Р'і Т К І Р' 1 +у-І)-а*(і,к 1 ) = Р'і К, д . 

Параметр регуляризации у выбирался путем ре- 
шения оптимизационной задачи 

У* = аг§тіп У (^НД-а^СуІ-аДсуІ-РДО )) 2 

7 

і-п-т 

методом дихотомии [8] (л=3). При п< 3 параметр у 
выбирался равным 0,01. 

Относительные ошибки оценок пластового да- 
вления и продуктивности скважины определялись 
по формулам: 

ад: (а = аъадо,) -аду аду 

У (« 2 * 0, Ж)) = аЬ§((а 2 *у , /?’) - а 2 ) / а 2 ), 1=1, п, 

где Р’ м (()=аІ(1,к")/а2(1,к'), а 2 ( 1 ,к') - оценки пласто- 
вого давления и продуктивности скважины, Р ы (і), 
а 2 - заданные (точные) значения пластового да- 
вления и продуктивности скважины (7). 

В качестве примера, в табл. 1, 2 приведены ре- 
зультаты расчета относительных ошибок оценок 
пластового давления и продуктивности скважины 
в зависимости от вида интегрированной модели (7) 
и объема измерений дебита скважины п. 

Для указания вида модели (7) введены следую- 
щие обозначения: {д‘} - априорная информация от- 
сутствует; {д\д}, {д",а}, {д\д,а}- учет промысловых 
данных о приемистости нагнетательной скважины, 
априорных данных и экспертных оценок продук- 
тивности скважины, промысловых данных о при- 
емистости и экспертных оценок продуктивности. 

Уровень ошибок дебита скважины составлял 
порядка 10 %, что соответствует точности наиболее 
распространенной измерительной системы типа 
«Спутник» [10]. Для имитации дебита скважины с 
относительным уровнем ошибок 10 % параметр с х 
выбирался равным 0,1. Для имитации дополни- 
тельных априорных сведений о пластовом давле- 
нии и продуктивности скважины с относительным 
уровнем ошибок 5 % параметры с 3 , с 4 в модели (7) 
выбирались равными 0,05. Данный уровень оши- 
бок соответствует приемлемому в практике нефте- 
газодобычи уровню погрешностей оценок гидро- 
динамических параметров нефтяных пластов [10]. 

Из табл. 1, 2 видно, что предложенные оценки 
параметров моделей дебита скважин с учетом до- 
полнительной информации о приемистости нагне- 
тательной скважины, дополнительных априорных 
сведений о пластовом давлении и продуктивности 
скважины (8) более точны по сравнению с оценка- 
ми метода наименьших квадратов (9). Значитель- 
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ный выигрыш в точности наблюдается при неболь- 
шом объеме данных истории разработки, порядка 
3-6 значений дебита и забойного давления сква- 
жины. Например, для получения оценок пластово- 
го давления с относительными ошибками в преде- 
лах 3 % достаточно пяти измерений дебита и забой- 
ного давления скважины (табл. 1). Для модели {д*}, 
не учитывающей дополнительные априорные дан- 
ные, относительная ошибка составляет 9 %. 


Таблица 1. Относительная ошибка оценок пластового давле- 
ния, % 


Модель 

дебита 

скважины 

Объем измерений дебита и забойного давления 
скважины 

і 

2 

3 

4 

5 

6 

7 

8 

9 

10 

{?} 

- 

49 

41 

21 

9 

7,1 

5,4 

4,7 

4,1 

3,3 

(<МІ 

41 

28 

18 

5,5 

5,3 

4,7 

4,1 

3,81 

3,4 

2,6 


39 

27 

16 

6,3 

4,8 

4,2 

3,6 

3,1 

2,9 

2,4 

{?’,?> а} 

27 

19 

5,7 

3,8 

2,9 

2,51 

2,34 

2,27 

2,16 

2,12 


В заключение отметим, что предложенные инте- 
грированные модели, метод и алгоритмы иденти- 
фикации дебита скважин (1), (4), (5), (8) позволяют: 
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Таблица 2. Относительная ошибка оценок продуктивности 
скважины, % 


Модель 

дебита 

скважины 

Объем измерений дебита и забойного давления 
скважины 

і 

2 

3 

4 

5 

6 

7 

8 

9 

10 

{</} 

- 

45 

32 

18 

11 

7,3 

5,3 

4,8 

4,6 

4,12 

<?,?} 

32 

23 

14 

8,3 

6,4 

5,9 

5,3 

4,7 

4,1 

3,6 

{4, а} 

28 

23 

9 

7,5 

5,3 

4,7 

4,3 

4,1 

3,6 

3,4 

а} 

21 

15 

4,9 

3,2 

2,9 

2,43 

2,12 

2,05 

2,01 

1,85 


1. Учитывать данные о дебитах (приемистости) 
скважин окружения, дополнительные априор- 
ные сведения и экспертные оценки гидродина- 
мических параметров (пластовое давление, 
продуктивность скважин и т. д.). 

2. Получать оценки в условиях априорной неопре- 
деленности о моделях зависимости дебита ис- 
следуемой скважины от дебита (приемистости) 
соседних скважин окружения, моделях допол- 
нительных априорных сведений и экспертных 
гидродинамических параметров. 

3. Значительно, более, чем в два раза, повысить 
точность оценок пластового давления и продук- 
тивности скважины. 
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